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ABSTRACT 

Historically,  peristalsis,  the  rhythmic,  traveling  wave  like 
pumping  of  bodily  organs  such  as  the  ureter,  has  been  mechanically 
simulated  by  the  extended  longitudinal  motion  of  rollers  squeezing 
fluid  along  a  flow  channel.   Blood  pumps  are  an  example.   The 
experimental  work  associated  with  this  thesis  has  developed  a 
peristaltic  pump  wherein  pumping  is  produced  by  the  limited 
transverse  motion  of  electromechanical  transducers.  The  transducers 
are  arranged  in  a  line  and  excitation  is  time  phased  so  that 
their  sequential  transverse  motions  produce  a  traveling  deflection 
wave  along  the  punp  flow  tube. 

The  pump  is  a  simple  device  made  from  transformer  cores. 
Each  transformer  serves  as  a  transducer;  the  "E",  with  windings 
on  the  center  leg,  provides  AC  and  DC  excitation;  the  "I"  serves 
as  the  actuator  which  squeezes  the  pump  tube  and  produces  flow. 
The  magnetic  action  is  much  like  a  solenoid  in  that  the  "E"  and 
"I"  serve  as  a  magnetic  circuit  and  the  magnetic  field  in  the  air 
gaps  produces  the  tube  squeezing  force. 

The  theory  developed  in  the  thesis  involves  the  electromechanical 
characteristics  of  each  pump  transducer  and  the  fluid  mechanical 
floitf  characteristics  for  the  pump.  The  final  tying  together  of 
theory  and  experiment  is  a  comparison  of  theoretical  pump  pressure- 
flow  characteristic  curves  and  experimental  pump  operating  points. 
Operating  points  do  conform  fairly  well  to  the  pressure-flow 
characteristics,  but  insufficient  operating  points  were  determined 
to  state  that  the  pressure-flow  curves  accurately  predict  pump 
operation. 

Production  of  peristaltic  pumping  by  transverse  deflections 
has  great  potential  value.   If  a  peristaltic  pump  is  ever  to 
be  used  as  a  replacement  for  an  organ  of  the  body,  longitudinal 
rollers  are  not  practical.  The  transverse  deflection  pump  is 
simple  and  also  might  be  miniaturized.  In  place  of  transformers, 
piezoelectric  crystals  could  be  used  as  transducers.  With  a 
DC  voltage  excitation,  crystals  could  produce  great  forces  over 
short  displacements.  Miniaturization  and  improvements  as  recommended 
in  the  thesis  are  worthwhile  areas  for  future  experimentation 
and  development. 
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SYMBOLS 

Notation  Name  and  Unit  of  the  Symbol 

a  Transducer  outer  leg  thickness,  meters 

b  Transducer  exciter  center  leg  thickness,  meters 

B  Magnetic  flux  density,  webers/m 

c  Phase  velocity,  m/sec 

d  Transducer  exciter  width,  meters 

D  Damping  coefficient 

V  Vector  differential  operator 

f  Frequency,  hertz 

f  Electric  force,  newtons 
e 

f  Fluid  pressure  force,  newtons 
P 

f  Spring  force,  newtons 
s 

ffi  Magnetic  flux,  webers 

p 

g  Gravity  acceleration,  m/sec 

h  Normalized  channel  width,  r/R 

H  Magnetic  field  intensity,  amp-turns /m 

H  Pump  suction  head,  meters 

H,  Pump  discharge  head,  meters 

i  Time  varying  current,  amps 

I   .  Steady  state  current,  amps 

3  (-1)* 

P 

J  Current  Density,  amps/m 

k  Wave  number ,  2  TT/  A 

K  Spring  constant,  newtons /m 
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Notation  Name  and  Unit  of  the  Symbol 

1  Actuator  width,  ty^,  meters 

L  Inductance,  henries 

X  Wavelength,  meters 

M  Actuator  mass,  kg 

jj.  Dynamic  viscosity,  ncwton-sec/m 

XCo  Permitivity  of  free  space,  henries/m 

M  Exciter  winding  turns,  turns 

V  Kinematic  viscosity,  M.J p 

P  Fluid  pressure,  newtons/m 


2 

m 


AP^  Pressure  rise  over  a  wavelength,  Mewtons/: 

Aj?  Normalized  pressure  rise,  ~£~r~jrz 

Q   ,;-:'      •  Instantaneous  fluid  flow,  m  /sec 

<QV  Time  average  fluid  flow,  m  /sec 

<Q>  Normalized  time  average  fluid  flow,  —= — 
c  Re 

r  Time  varying  component  of  channel  width,  meters 

R  Equilibrium  channel  width,  meters 

P  Fluid  density,  kg/m 

t  Time,  sec 


T  Maxwell  stress  tensor,  Newtons/m' 

v  Fluid  velocity,  m/sec 

w  Channel  width,  meters 

W  System  energy,  watts 

W  System  co-energy,  watts 

x  Spatial  co-ordinate 

7. 
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Notation  Name  and  Unit  of  the  Symbol 

X  Transducer  gap  width,  meters 

X  Tube  support  thickness,  meters 

y  Spatial  co-ordinate 

z  Spatial  co-ordinate 

CJ  Frequency,  2Tff,  radians/sec 
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CHAPTER  I 
INTRODUCTION 

Peristalsis  is  the  rhythmic  motion  of  the  walls  of  the 
alimentary  canal,  the  ureter,  and  other  hollow  organs  of  the 
body,  consisting  of  traveling  wave  contractions  and  dilations 
that  move  the  contents  of  the  tube  onward.   Outside  the  body, 
peristaltic  pumping  can  be  achieved  by  mechanically  squeezing 
a  fluid  filled  channel  or  tube  in  a  traveling  wave  manner. 

The  potential  value  of  such  a  devi.ce  must  be  apparent. 
In  industry,  the  need  often  exists  to  pump  fluids  without 
having  the  fluid  touch  the  pump  apparatus,  as  in  the  case 
of  blood  and  caustic  chemicals.  The  peristaltic  pump  enables 
the  insertion  of  a  fluid  carrying  tube  into  the  pump  with 
no  actual  contact  between  the  fluid  and  the  pump  mechanism. 
In  addition,  if  such  a  device  could  be  miniaturized,  it  might 
be  inserted  into  the  body  to  create  artificial  body  functions 
when  natural  functions  have  failed. 

All  present  peristaltic  pumps  rely  on  transverse  pump 

tube  deflections,  which  travel  longitudinally,  to  create  fluid 

(12  3) 

motion.        Deflection  is  generally  produced  by  a  mech- 
anical roller  which  moves  along  the  tube,  squeezes  it,  and 
pushes  the  fluid  onward.   Electrical  forces  are  involved  only 
in  that  an  electric  motor  produces  the  force  which  causes 
the  longitudinal  mechanical  motion. 
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This  investigation  deals  with  a  different  class  of  electro- 
mechanical devices.  In  general,  mechanical  motion  is  purely 
transverse  and  electric  forces  are  directly  related  to  pump 
tube  deflections.     Schematic  representations  of  such  devices 
are  shown  in  Figure  1„1.  In  Figure  1.1. a,  the  channel  walls 
form  a  continuum  electromechanical  transducer.  The  source 
of  electrical  force  is  a  current  sheet  which  is  wound  similar 
to  the  windings  in  an  electric  motor  field  and  which  produces 
a  sinusoidal  traveling  wave.  If  the  channel  is  an  elastic  tube, 
the  continuum  transducer  can  be  approximated  by  arranging 
a  large  number  of  discrete  magnetic  segments  on  top  of  the 
tube,  Their  response  to  the  traveling  magnetic  wave  will 
cause  tube  deflections. 

Figure  1,1, b  represents  a  peristaltic  pump  which  also 
moves  purely  transversely,  but  the  sources  of  tube  deflection 
are  discrete  transducers  rather  than  an  electromechanical 
continuum.   The  schematic  shows  each  transducer  being  one 
quarter  wavelength.  In  order  to  produce  a  single  wavelength 
of  a  traveling  wave,  as  illustrated,  each  transducer  excitation 
must  be  phase  shifted  ninety  degrees  from  its  predecessor. 
The  number  of  discrete  transducers  per  wavelength  can  be  in- 
creased as  desired,  provided  they  are  properly  excited.  In 
the  limit  of  many  transducers,  the  discrete  pump  approaches 
being  a  continuum. 


10. 


SCHEMATIC  REPRESENTATION  OF  A  TRAVELING  WAVE  PERISTALTIC  PUMP 


Current  Sheer 


a. 


A 


Transducer 
Actuator  % 


JZ3IL 


xz^zi 


b. 


Transducer  Exciters  « 

A 


_V^sa» 


A» 


R 

Pressure  =  P(z) 
Velocity  -  V  (x) 
r  -   r(z) 


Figure  1 . 1 
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In  both  representations,  Figures  1.1. a  and  1,1, b,  the 
electric  force  which  causes  tube  deflection  is  a  function 
of  the  gap  width,  X,  between  the  actuator  and  the  exciter. 
It  will  be  shown  that  this  relationship  produces  a  mechanical 
instability  which  causes  the  discrete  configuration  of  Figure 
ld.b  to  operate  in  an  on -off,  that  is,  channel  open-channel 
cut  off,  mode. 

The  analytical  model  for  the  pump  is  based  upon  the  physical 

situation  most  likely  to  be  encountered  where  peristaltic 

(5) 

flow  exists.     Transverse  tube  dimensions  are  small  compared 

to  a  wavelength  of  the  traveling  wave0  Transverse  velocities 
and  pressure  gradients  are  small  compared  with  longitudinal 
ones.  The  pressure  in  the  tube  is  assumed  to  be  instantaneously 
uniform  over  each  cross  section0  The  fluid  density  and 
viscosity  are  assumed  to  be  constant.  The  Reynold's  number 
is  assumed  to  be  low  enough  so  that  flow  is  viscous.  In  order 
to  insure  the  validity  of  a  viscous  model,  Reynold's  numbers 
should  be  less  than  one0 

The  tube  walls  move  in  the  transverse  direct] on  only. 

The  traveling  wave  shape  will  be  analyzed  as  both  sinusoidal 

(7  8^ 

ands.quare,   '    Pumping  does  not  seem  to  be  sensitive  to  the 

shape  of  the  wave,  but  specific  values  of  flow  rates  are  sensitive 
to  the  type  of  tube  deflection. 

The  pump  developed  for  investigation  in  this  work  is 
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shown  in  Figure  1.2.  It  consists  of  four  transducers  which 
are  excited  to  produce  a  single  traveling  wavelength  along  the 
pump  tube.  The  tube  is  inserted  between"  the  transducer  exciter 
and  actuator.  The  actuator  is  supported  by  four  springs, 
one  at  each  corner,  and  by  the  fluid  pressure  in  the  tube 
passing  underneath.  The  transducers  have  two  exciter  windings. 
One  supplies  direct  bias  cut-rent.  The  other  supplies  alter- 
nating current.  The  alternating  current  source  is  a  sine- 
cosine  potentiometer,  the  output  of  which  is  amplified  through 
amplidynes  to  produce  the  desired  lox^  frequency  current  values. 
The  phase  shift  in  excitation  between  each  transducer  is  ninety 
degrees. 

The  analysis  of  the  pump  is  done  step  by  step,  starting 
with  a  single  transducer  and  ending  with  the  completed  per- 
istaltic pump. 

In  chapter  two,  the  inductance  of  a  single  transducer, 
the  supporting  spring  constant,  and  the  electric  forces  exerted 
on  a  transducer  as  it  will  appear  in  the  pump  are  determined. 
Then  the  time  variation  of  a  single  transducer  and  the  set 
of  four  transducers  in  sequence  is  ascertained  for  various 
direct  and  alternating  current  excitations.  The  frequency 
of  alternating  current  is  arbitrary  except  that  the  Reynold's 
number  is  directly  proportional  to  frequency.  Lower  Reynold's 
numbers  result  from  lower  frequencies.  The  lowest  possible 
Reynold's  number  is  desired  to  have  viscous  flow  forces  dominate 
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PERISTALTIC  PUMP  PARAMETERS 


a  «=  1.45  X  10~2meters 

b  «=  2o90  X  10_2rneters 

d  -  2.95  X  I0~2meters 

1  «=  7oOO  X  10~2meters 


N  .     =  500  turns 
primary 

N     ,    =26  turns 
secondary 

Xs  -  U73  X  10"\eters 
w  «=  variable  tube  width 
X  =  variable  gap  width 


Figure  1.3 
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inertial  flow  forces. 

Once  the  transducers  are  analyzed  individually,  the  per- 
istaltic pump  of  which  they  are  a  part  is  analyzed  as  a  whole 
in  chapter  three 0  First,  an  empirical  expression  for  the 
pressure  force  exerted  by  the  pump  tube  is  obtained.  Then 
a  theoretical  expression  describing  transducer  actuator  motion 
is  determined.  Finally,  based  on  known  actuator  deflections, 
the  pressure-flow  characteristics  of  the  pump  are  determined 
for  a  square  traveling  wave  and  a  sinusoidal  traveling  wave 
description  of  tube  deflections.  The  experimental  operating 
points  for  the  pump  are  compared  to  these  characteristics. 

Throughout  chapters  two  and  three,  analytic  work  is 
immediately  followed  by  test  results.  This  is  done  because 
experimental  results  from  some  sections  are  used  in  following 
analytic  work. 
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CHAPTER  II 
MAGNETIC  TRANSDUCER  ANALYSIS 
2,1  Magnetic  Transducer  Inductance 

2.1.1  Theoretical  Inductance 

The  theoretical  inductance  of  a  single  transducer  is  det- 

(9) 

ermined  us:mg  lumped  parameter  electromechanical  techniques. 

Apply  Ampere's  law  to  the  transducer  as  pictured  in  Figure  2.1a. 
The  permeability  of  the  exciter  and  actuator  are  assumed  to  be 
infinite.  This  requires  that  the  magnetic  field  intensity,  H, 
must  be  zero  since  the  flux  density,  B,  is  finite  and  B  =/<_H, 
The  only  magnetic  field  Is  in  the  gaps  between  the  exciter  and 
the  actuator, 

j  7.XH  ■  nda=  $s  J  •  n  da  2.1a 

Apply  Stokes*    theorem  to  obtain 

j ,  H   •   dl  -    \    J    •   n  da  2.1b 

The  contours  for  these  integrals  are  indicated  as  dotted 
lines  on  Figure  2,1a.   Integration  yields 


Hj  X  +  H£  X  =  N  I  2.2 


Taking  a  surface  that  encloses  the  transducer  actuator  in 
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Figure  2.1 
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Figure  2,1b,  knowing  that  V*  B  --=  0,  integrating  V*  B  over  the 
enclosed  volume,  and  applying  the  divergence  theorem  yields 

\,  B  •  n  da  =  0  2.3 

In  terms  of  the  magnetic  field  penetrating  the  enclosing  surface, 
this  becomes 

2ij)H1ad-jU[)H  bd=0  2.4 

Solve  equations  2.2  and  2,4  simultaneously. 


H   _  _b N_I  , 

hl   "   2a  +  b  X  *° 


H 2a__  N_I  9  , 

H2  ~  2a  +  b  X  2'6 


The  total  flux  through  the  center  leg  is 


j>  =  MQ   H2  b  d  2.7 


2^0  N  1   a  b  d  2.8 

X  (2a  +  b) 


In  the  absence  of  leakage,  the  flux  links  the  N  turn  winding 
N  times.   The  flux  linkage  is 
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X  =  N  <t>  2.9 


2.10 


2  XL  1 


X  (2a  +  b) 
=  L(X)  I  2.11 


L(X)  -  2^0  N"  a  b  d  2.12 

X  (2a  +  b) 


2.1.2  Experimentally  Measured  Inductance 
The  transducer  exciter  has  a  primary  winding  of  500  turns 
and  a  secondary  winding  of  26  turns  on  the  center  leg,  A  com- 
parison is  to  be  made  between  the  theoretical  transducer  inductance 
and  the  measured  primary  and  secondary  inductances.   Because  the 
secondary  inductance  winding  has  only  26  turns,  it  must  be  nor- 
malized in  proportion  to  the  turns  ratio  squared  if  it  is  to  be 
made  comparable  to  the  primary  inductance  on  the  same  plot  of 
inductance  versus  gap  x^idth, 

L(X)       -,.   ,  =  (  ^  )2  L(X)         .       2.13a 
x  's  normalized      26        s  measured 


=  370  L(X)  2.13b 

sm 


The  theoretical  transducer  inductance  is  evaluated  from 
equation  2.12  by  substituting  the  transducer  parameters  from 
Figure  1.3. 
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Plots  of  transducer  inductances  versus  gap  width  are  shown 
in  Figures  2.3  and  2,h,     Further  comparative  insight  can  be 
gained  by  taking  the  natural,  logarithm  of  equation  2.14  and 
differentiating  with  respect  to  In  X, 


In  L(X)t  *=  In  1.34  x  10"^  -  In  X  2.15a 


d  ln  h™t   .  -1  2.l5b 

d  In  X 


A  logarithmic  plot  of  theoretical  transducer  inductance 
versus  gap  width  is  a  straight  line  with  slope  of  minus  one. 
Figure  2,5  shows  logarithmic  plots  of  theoretical,  primary, 
and  normalized  secondary  inductances. 

The  measured  values  of  primary  and  secondary  inductance  vary 
considerably  from  the  theoretical  values.  Inspection  of  Figure 
2,11  shows  that  electric  forces  calculated  using  the  theoretical 
transducer  inductance  are  closer  to  the  actual  value  of  measured 
electric  force  than  electric  forces  calculated  using  measured 
inductance  values.  The  range  of  gap  width  for  these  calculations 
is  from  .003  to  .00^-  meters.  In  the  limiting  cases  of  X  =  0  and 
X  -  infinity, the  measured  inductance  values  are  more  accurate 
than  theoretical  ones.  Theoretically,  as  X  approaches  infinity, 
the  inductance  approaches  zero.  And,  as  X  approaches  zero,  the 
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EXPANDED  PLOT  OF  PRIMARY  AND  THEORETICAL'  TRANSDUCER  INDUCTANCES 
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inductance  becomes  infinite.   Practically,  neither  of  these  the- 
oretical limits  is  accurate.  In  the  real  infinite  limit,  the 
gap  width  is  so  great  that  the  magnetic  field  lines  find  a  path 
of  less  reluctance  by  going  directly  from  the  center  exciter  leg 
to  the  outside  legs  through  the  air.  Even  in  the  case  of  less 
than  infinite  gap  width,  some  field  lines  travel  directly  from 
leg  to  leg  altering  the  transducer  inductance.   In  the  zero  limit, 
when  no  air  gap  exists  between  the  actuator  and  the  exciter,  the 
actual  inductance  is  not  infinite  because  the  permeability  of 
the  actuator  and  exciter  metal  is  not  infinite,  as  was  assumed 
in  the  theoretical  model,  but  finite. 

Fringing  fields  and  non-infinite  transducer  permeabilities 
affect  the  measured  inductance  values  for  all  gap  widths.  Field 
losses  can  account  for  the  difference  between  primary  and  sec- 
ondary inductances.   The  primary  exciter  winding  is  wound,  directly 
to  the  exciter  center  leg.   The  secondary  winding  is  wound  over  the 
primary.   The  tendency  would  be  for  greater  losses  to  be  experienced 
in  the  secondary,  resulting  in  a  larger  inductance.  Figure  2.3 
shows  this  to  be  true. 

Additional  errors  could  have  come  from  the  impedamce  bridge 
used  for  inductance  measurements, 
2.2  Actuator  Support  Spring  Constant  ■  . 

Each  transducer  actuator  is  supported  by  four  springs.  The 
value  of  the  spring  constant,  K,  for  these  springs  was  desired. 
Weights  were  placed  on  top  of  each  actuator.  The  resulting 
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PLOT  OF  SPRING  FORCE  ACTING  ON  THE  TRANSDUCER  ACTUATOR  VERSUS  SPRING 

COMPRESS  I  UN 
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Figure  2.? 
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change  in  gap  width  due  to  these  weights,  which  is  the  spring 
deflection,  was  measured  with  micrometers. 

The  values  of  weights  and  spring  deflections  are  plotted 
on  Figure  2.7.  The  spring  constant  for  the  four  springs  acting 
in  parallel  is  the  slope  of  the  curve.  Thus,  the  spring  force, 
f  ,  is  expressed  as  K  (Xn  -  X)  where  X  is  the  undeflected  gap 
width.  Figure  2.7  shows  the  spring  constant  to  be  2860  newtons/ 
meter, 
2.3  Electric  Forces 

2.3*1  Electric  Force  Determination  Using  Energy  Methods 
The  magnetic  transducer  is  represented  as  a  two  terminal 
circuit  in  Figure  2.8a.  The  total  energy  stored  in  the  system 
is  W  .  The  power  or  rate  of  change  of  energy  is 


1  =  IV  -  I   TT  2. ID 

e  dt 


dt 


=  i^f  -  f  %  2.17 

dt     e  dt 


dW  =idX-f  dx  2.18 

m  e 


It  is  desired  to  have  current  as  the  variable,  so  a  co- 
energy  is  defined  as 


d  W  =  -d  W  +  i  X  2.19 

m       m 
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Substitute  equation  .2.19  in  2,l8. 


d  W  -  \d±   +  f  dx  2.20 

m   '       e 


Integrate  the  co-energy  along  the  path  shown  nri  Figure  2.8b, 
The  integral  is  independent  of  the  path  taken.  Along  path  (Ij  , 
current  is  zero  so  the  electric  force  is  zero.  Along  path  (?)  , 
the  displacement  is  at  its  final  value,  X. 

W  =  \*f  dx  +  51  X  di  2.21 

m   Jo  e      -'o 

c--o     X  =  X 

Using  the  value  of  X  from  equation  2.10  and  the  definition 
of  L(X)  from  2.12, 


W  -  3   N^_a_b  d  I2  2.22 

m     X~(2a  +  b)~ 


|  L(X)  I2  2.23 


In  differential  form,  d  W  can  be  represented  as 

m 


d  W 


•J*  = m  di  +  m  dx  2.2'+ 


di        dx 


Compare  this  equation  with  equation  2.20, 


di 

31. 


*B-iS  2.26 

6   dx 


From  equation  2.22,  the  electric  force  is 


f  =.4N2abd     2 
e    _ 

T'   (2a  +  b) 

_  i  d_l(X)   T2 
~  2  dX   •  J- 


2.27 


2.28 


2»3»2  Electric  Force  Determination  Using  Maxwell 

04.         rn         (lO 

Stress  Tensor 
The  electric  force  on  the  transducer  actuator  in  the  x  direction 
is  desired.  This  is  evaluated  using  the  definition  of  the  stress 
tensor  in  Figure  2.9.  Since  the  only  magnetic  fields  are  in  the 
x  direction,  only  an  x  component  of  the  stress  tensor  exists. 

Txx*=*VHl+H2  +  Hl  >  2-29 

fe  "  k  TxX  '  »  da  2-30 

■=  i  -M0  (  2  H2  a  d  +  H2  b  d  )  2.31 

Substitute  equations  2.5  and  2.6  for  H  and  H„. 

f  =  ^0  ^  *  b  d  T2  2.32 

e      ! 

X^  (2a  +  b) 

which  is  identical  to  equation  2,2?. 
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MAXWELL  STRESS  TENSOR  REPRESENTATION  OF  THE  TRANSDUCER  ACTUATOR 

ELECTRIC  FORCE 
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2.3 t 3  Experimental  Evaluation  Of  Electric  Force 
For  any  applied  current  to  the  excite.r  windings,  the  actuator 
will  assume  an  equilibrium  position  such  that  the  springs  are 
compressed  and  exert  an  upward  spring  force  which  is  balanced 

by  the  applied  electric  force  and  the  weight  of  the  act\iator 

4  •    ,      j  (12) 
acting  d  ownwa r d , 

The  actuator  weight  is  balanced  by  an  equilibrium  spring 

f orce0 

Kg  =  K  (  XQ  -  X1  )  2o33 

The  electric  force  will  then  be  balanced  by  the  additional 
compression  of  the  springs. 

f   .  -  K  (  X.  -  X  )  2.3^ 

era!        1 

This  force  is  the  experimentally  measured  electric  force. 

The  electric  force  can  be  determined  in  two  other  ways. 
Theoretically  it  is 


i  d  L(X) 
etl 


f  -  =  *    •   t 


dX 


l2  2.35 


X 


from  equation  2.28  where  ,v  *    't  is  the  slope  of  Figure  2,/4-. 

x  ■  OA 

By  use  of  measured  inductances  and  the  slope  of  the  measured 
inductance  curve,  measured  electric  forces  are 
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PLOT  OF  SPRING,  THEORETICAL  ELECTRIC,  AND  MEASURED  ELECTRIC  FORCES 

VERSUS  GAP  WIDTH 
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Over  the  range  of  gap  widths  measured,  the  theoretical  electric 
force  most  accurately  predicts  the  actual  electric  force,  f  , . 
The  electric  force  calculated  from  the  measured  inductance  is 
considerably  in  error.  It  was  evaluated  using  the  slope  of  the 
measured  inductance  curve,  Figure  2o^-,  from  section  2,1,2.  An 
error  might  have  been  expected  for  the  reasons  previously 
mentioned  in  section  2,1,2, 
2.^-  Magnetic  Transducer  Force  Equilibrium 

Newton's  law  applied  to  the  transducer  actuator  is 


A2  A 

M  —  +  D  —  =  f   .   +  f  -,   .  .   -  Mg  2.37 

,,2     ,,     spring    electric 


The  position  that  the  transducer  actuator  takes  in  steady 

state  equilibrium  represents  a  balance  between  the  upward  spring 

(13) 
force  and  the  downward  electric  and  weight  forces. 


F(X)  =  f  +  f  -  Mg  2.38 


2 

=  K  (  X0  -  X  )  +  ^0  li     a  b  d  I2  -  Mg        2.39 

X2  (2a  +  b) 


2      2 

K  (Xq-X^  +  K  (Xr  X)  +  M0   M  abd  J~-  Mg    2.4 


X2  (2a+b) 


The  weight  of  the  actuator  cancels  K  (XQ  -  X  ), 

37. 


F(X)  -  K  (Xx  -  X)  4  -^0  N  *  abd  *  2.«H 

X2  (2a  +  b) 

Figure  2.13  is  a  plot  of  these  two  functions  of  X.  The 
points  where  they  intersect  are  the  equilibrium  operating  points 
or  the  equilibrium  positions  assumed  by  the  actuator.  A  separate 
plot  will  be  obtained  for  any  value  of  current.  The  points,  X  , 
are  unstable  equilibrium  points,  A  change  in  X  from  this  value 
causes  a  force  imbalance  which  tends  to  cause  X  to  change  even 
further.  The  points,  X  ,  are  stable.  Any  change  in  X  from  X 
causes  a  force  imbalance  which  tends  to  return  the  actuator  back 
to  X  o  The  actuator  always  takes  a  stable  position. 

As  current  is  increased  beyond  0,677  amperes,  the  values 
of  electric  force  are  always  greater  than  spring  force.  There 
are  no  equilibrium  positions  for  the  actuator.  Electric  force 
always  overcomes  the  springs  and  closes  off  the  gap.  The  point 
of  tangency,  Xrp,  is  the  equilibrium  position  for  the  actuator 
when  current  is  equal  to  0.677  amperes.  Any  perturbation  acting 
to  increase  X  causes  the  actuator  to  return  to  X  .   But,  any 
perturbation  which  acts  to  decrease  X  will  cause  the  actuator 
to  close  off  to  X  =0. 

The  accuracy  of  Figure  2.13  can  be  checked.  Each  stable 
equilibrium  point,  X  ,  on  Figure  2.13  should  correspond  to  the 
measured  gap  width  from  figure  2,10  for  the  same  value  of  current. 
Examination  shows  this  to  bo  true. 
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PLOT  OF  SPRING  AMD  ELECTRIC  FORCES  VERSUS  GAP  WIDTH  FOR  VARIOUS  CURRENTS 
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Figure  2.13 
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2.5  Time  Variation  Of  Transducer  Actuator  Position  For  Sinusoidal 
Applied  Current 

Plots  of  actuator  position  versus  time  are  obtained  using 
the  current-gap  width  correspondence  from  Figure  2.10.  Figures 
2.1^a,  2.15a,  2.l6a,  and  2,17a  are  plots  of  time  varying  applied 
current  and  time  varying  head  position  for  various  combinations 
of  direct  bias  current  and  alternating  current.  The  gap  width 
"closes  off"  at  X  -  .00173  meters  because  the  tube  support  spacer, 
as  shown  in  Figure  1.2c,  prevents  it  from  going  to  zero.  When 
the  pump  tube  is  inserted,  it  will  be  squeeze'd  completely  shut 
in  this  condition. 

It  was  desired  to  operate  the  pump  less  than  fully  occluded. 
In  concept  this  could  have  been  accomplished  by  limiting  the  exciter 
currents  such  that  the  actuator  never  fully  closed  off  the  channel. 
Such  operation  is  illustrated  in  Figure  2,17a.   Under  these  conditions, 
the  channel  has  a  steady  state  width  and  a  time  variation  in 
width  around  the  steady  state.  Schematically  this  is  shown  in 
Figure  1,1  when  R>r,  Practically,  operation  of  the  pump  as 
described  is  very  difficult.  During  operation,  variation  in 
the  characteristics  of  the  individual  transducers  from  each 
other  cause  one  to  operate  properly  while  another  closes  off. 
Figure  2,17c  shows  actuator  three  operating  almost  properly, 
while  the  other  three  are  closing  off.  The  transducers  are  so 
unstable  that  only  slight  perturbations  are  necessary  to  make 
the  gap  close  off  completely. 

fcl. 


To  overcome  this  difficulty,  a  channel  limiter  was  placed 
astride  the  pump  tube.   Actuator  motion  is  restricted  because  the 
actuator  contacts  the  rigid  limiter  as  it  tries  to  c3 ose  off. 
The  pump  transducers  can  be  driven  as  hard  as  desired,  but  peristalsis 
takes  place  only  on  that  section  of  the  tube  above  the  limiter. 
The  thickness  of  the  limiter  can  be  varied  to  obtain  the  desired 
values  of  R  and  r.   Figure  2.1 6a  illustrates  the  time  variation 
of  actuator  motion  with  the  limiter  in  place  and  Figure  2,l6c 
shows  the  Sanborn  recordings  of  actuator  motion  under  these 
circumstances , 

Figures  2.1^4-b,  2.15b,  2.16b,  and  2.17b  are  representations 
of  the  positions  of  the  four  actuators  for  various  times  covering 
one  complete  cycle.   Each  actuator  is  phased  ninety  degrees  from 
the  preceeding  one0  The  result  is  a  traveling  wave  to  the  right 
which  can  be  approximated  as  a  square  wave  or  as  sinusoidal. 

Finally,  Figures  2.l4c,  2.15c,  2.1 6c,  and  2.17c,  two  of 
which  have  already  been  discussed,  are  Sanborn  recorder  plots 
of  actuator  position  versus  'time  for  all  four  transducer  actuators. 
The  plots  are  from  the  output  of  photo  electric  cells  which  were 
arranged  to  sense  the  motions  of  the  actuators.  Determination 
of  the  exact  magnitudes  of  gap  widths  is  not  possible  from  the 
recordings  because  the  recording  amplitudes  are  relative  and 
a  function  of  the  sensitivity  setting.  However,  the  recordings 
do  illustrate  the  phase  shift  between  the  actuators  which  creates 
the  traveling  wave.  The  phase  shift  is  not  ninety  degrees  between 
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each  actuator.   The  outputs  of  the  power  supplies  for  the  exciter 
currents  were  not  quite  properly  phased.  The  resulting  traveling 
wave  was  uneven:  close  close  pause  close  close  pause  close 
close. • . . . o • • • o  o 

A  comparison  betx^een  the  predicted  actuator  time  varying 
position  (e.g.  Figure  2.15a)  and  one  cycle  of  an  actual  measared 
time  variation  (e.g.  Figure  2.15c)  shows  that  the  actual  motion 
conforms  quite  well  to  predicted  motion.  The  operation  of  the 
transducer  is  effectively  open-closed.  The  transducer  is  so 
unstable  that  it  acts  almost  like  an  on-cff  relay0  This  leads  to 
the  conclusion  that  an  elaborate  multi-phase  sinusoidal  power 
supply  is  not  necessary  for  efficient  operation  of  the  pump0 
A  simple  on-off  direct  current  power  supply  for  each  transducer 
which  is  actuated  ninety  degrees  ahead  of  its  predecessor  is 
sufficient o 
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THEORETICAL  TRANSDUCER  ACTUATOR  POSITIONS  FOR  VARIOUS  TIMES 
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CHAPTER  III 
PERISTALTIC  PUMP  ANALYSIS 
Analysis  of  the  electrical  and  mechanical  properties  of 
individual  transducers  has  been  completed  in  the  previous  chapter.  ' 
This  chapter  is  concerned  with  the  analysis  of  the  peristaltic 
pump  as  a  whole.   The  fluid  filled  tube  running  between  each 
transducer  actuator  and  exciter  exerts  an  upward  pressure  force 
on  the  actuator  that  must  be  determined  and  entered  into  the 
force  balance  on  the  actuator.  The  magnitude  of  this  upward 
pressure  force  relative  to  the  electric  force,  which  acts  to 
pull  the  actuator  downward  and  close  off  the  flow  channel,  must 
also  be  determined  for  the  following  reasons.   If  the  pressure 
force  is  comparable  to  or  greater  than  the  electric  force,  the 
periodic  motions  of  the  transducer  actuators  (i.e.  tube  deflections) 
are  not  straight  forward  to  analyze.  Fluid  mechanical  effects 
between  transducers  must  be  accounted  for.   If  the  electric  force 
dominates,  the  tube  deflections  have  the  same  periodicity  as  the 
applied  magnetic  traveling  wave  and  approximate  analysis  of  tube 
deflections  is  simple. 
3..1  Pump  Fluid  Pressure  Force 

3*1.1  Theoretical  Description  Of  The  Pressure  Force 
Exerted  On  The  Transducer  Actuator  By  The  Pump 
Suction.  Head 
The  static  fluid  pressure  exerted  on  the  tube  wall  and  the 
transducer  actuator  is  a  function  of  the  suction  head,  H   Inertial 
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or  fluid  velocity  effects  on  the  fluid  pressure  will  be  neglected 
in  the  analysis. 

P-/?gHs  3.1 

The  total  fluid  pressure  force  acting  on  each  actuator  is 


f  -Pwl  3o2 


-/gwlHfi  3.3 

where  w  and  1  are  dimensions  from  Figure  1.2. 
Solve  equation  3. 3  £°?   wl. 

wl  -  J^  3.4 

wl  is  the  variable  area  over  which  the  fluid  pressure  acts 
on  the  actuator.  It  is  a  function  of  the  gap  x^idth  and  will 
be  derived  experimentally  in  the  next  section.  The  emperical 
expression  for  wl  would  be  used  in  describing  the  motion  of  the 
peristaltic  pump,  section  3.2,  under  operating  conditions  where 
the  fluid  pressure  is  comparable  to  the  electric  force. 

3.1.2  Experimental  Evaluation  Of  The  Fluid  Pressure 

Force 
Figure  3.2  is  a  plot  of  the  total  force  due  to  springs  and 
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PLOT  OF  SPRING  FORCE,   FLUID  PRESSURE  FORCE,    AND  COMBINED  SPRING-FLUID 
PRESSURE  FORCE  VERSUS  GAP  WIDTH 
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fluid  pressure  acting  on  each  actuator  versus  channel  width. 
Since  the  spring  force  is  known,  it  is  subtracted  from  the  total 
force,  leaving  the  force  due  to  the  fluid  pressure  alone.  The 
fluid  pressure  curves  yield  an  empirical  expression  for  the  fluid 
pressure  force  as  a  function  of  suction  head  and  channel  width. 


fp  =  1.2  +  6  Hg  -  ?50  (X  -  Xg)  3.5 


where  H  and  X  are  in  meters  and  f  is  newtons, 
s  p 

Substitution  of  this  empirical  expression  in  equation  3.^ 
yields 


1.2  +  6  H  -  750  (X  -  X  ) 

wl  = PSX "  3-6 


3.1.3  Comparison  Of  Relative  Magnitudes  Of  Forces  Acting 
On  The  Transducer  Actuator . 
The  forces  acting  on  the  actuator  vary  with  gap  width,  suction 
head,  and  exciter  current.  The  fluid  and  spring  forces  act  upward, 
while  the  weight  of  the  actuator  and  the  electric  force  act 
downward.  During  operation  the  gap  width  varies  from  about  ,00^ 
to  .0017  meters.  Table  2>.\    shows  values  of  forces  for  different 
gap  widths,  suction  heads,  and  exciter  currents. 

The  figures  in  Table  3*1  do  not  clearly  indicate  what  happens 
when  the  pump  operates.  As  current  to  the  exciter  reaches  its 
peak  value,  the  actuator  begins  to  close  due  to  a  net  imbalance 
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Table  3.1 
Approximate  Force  Values  For  Various  Gap  Widths 


Force, (newtons) 

X  =   .004 

X  =  .003 

X  =  .002 

X   =   .0017 

Fluid  Pressure  Force 
f    ,   Figure  3.2 

H  -0.1 
s 

H  =0.-2 

V0-3 

H  =0.4 
s 

0.10 
0.65 
1.25 
1.85 

0.80 
1.35 
1.95 
2.55 

1.50 
2.05 
2.70 
3.30 

1.70 
2.25 
2.90 
3.50 

Electric  Force 
f   ,   Figure  2.11 

1=1.02 

1=1.50 

0.42 
0.91 

3.60 
7.75 

8.05 
17.50 

11.10 
24.00 

Spring  Force 
f   ,   Figure  2.7 

2.80 

5.70 

8.60 

9.40 

Actuator  Weight 
Mg 

2.91 

2.91 

2.91 

2.91 

of  forces  acting  downward.  As  the  actuator  closes,  the  electric 

force  increases  at  a  rate  greater  than  the  spring  and  pressure 

forces  combined.   The  actuator  is  unstable.  Once  it  starts  to 

close,  the  imbalance  of  forces  increases  and  it  will  close  completely. 

Further,  as  the  channel  pinches  closed,  all  fluid  is  forced  out 

and  the  fluid  pressure  drops  to  zero.   Observation  of  operation 

of  the  pump  shows  that  the  fluid  pressure  and  fluid  dynamics 

do  effect  the  pump,  but  that  the  electric  force  does  dominate 

and  the  periodic  motion  of  the  actuators  has  the  same  frequency 

as  the  applied  exciter  current. 

(l4  15) 
3.2  Theoretical  Description  Of  Peristaltic  Pump  Motion   ' 

The  balance  of  forces  acting  on  the  actuators  in  Figure  1.2 
yields  6l . 


P  (wl)  +  K  (XQ  -  X)  =  Kg  +  f^  3.7 


o 


In  static  equilibrium,  the  weight  is  balanced  by  an  equilibrium 
spring  force. 


K  (XQ  -  X^  -  Mg  3.8 


Thus 


P  (wl)  +  K  (X1  -  X)  =  f e  3.9 


Rearrange  equation  3°9  and  divide  by  wl. 


p  =  £©  _  K  (X1  -  X) 

wl   z 

wl 


3.io 


Substitute  for  f  from  equation  2,2?  and  for  wl  from  equation  3«6. 


e 


Pg   Hs  f      2    2  T2  1 

p  -  urr^r^m  [  A> N  abd  J2  - K  (xi  -  X)J  3-u 

Consider  the  tube  wall  to  be  a  continuously  deforming  surface 
excited  by  a  traveling  wave.  That  is,  the  discrete  motions  of  the 
transducer  actuators  are  approximated  as  an  appropriate  continuous 
traveling  wave. 

The  tube  wall  deflection  travels  to  the  right  in  Figure  3.3 

with  a  constant  velocity,  c,  which  is  a  function  of  the  frequency 

and  wavelength  of  the  traveling  wave, 
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CROSS  SECTIONAL  VIEW  OF  THE  PERISTALTIC  PUMP  IN  THE  FIXED  AND 

MOVING  FRAMES 
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Figure  3^3 
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c  -  >f  3.12 

where  ^  =  41  frorc  Figure  1.2 

The  tube  wall  itself  is  fixed  in  longitudinal  position  and 
moves  only  transversely.  The  tube  wall  position  in  the  wave  or 
moving  frame  of  reference  is 

X  =  R  +  r(z')  3.13 

In  the  wave  frame,  the  traveling  wave  is  a  periodic  stationary 
deflection.   By  assuming  steady  state  flow  conditions  and  ignoring 
inertia,  the  Navier-Stokes  fluid  flow  equation  is 


Vp  =  jjisj2  y  3.1^ 


The  model  for  the  tube  or  channel  is  quasi  one  dimensional. 
The  pressure  is  independent  of  transverse  position.  The  velocity 
is  purely  in  the  longitudinal  direction,  but  varies  with  transverse 
position. 


d  P' 

S-i-  -  0  3.15 

dx' 

d  P'      d2  v'  _  ., 

z  3.1° 


dz'    M 


dx' 


In  the  wave  frame  of  reference,  the  volume  rate  of  flow 
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per  unit  depth  in  the  y  direction  is 

ft 

Q*  =  J   v*  dx  3,17  - 

Jo   z 

The  fluid  velocity  is  determined  by  integrating  equation 
3.l6  and  applying  the  boundary  conditions 


v»    (x  =  X)    =  vw   (x  -  X)    -  c  =  -c  3.18 


v«    (x  =  0)    =  -c  3.19 


The  result  is 


K  -  tic  f7  (x" " Xx)  " c  3-20 


Substitute  equation  3 .20  in  equation  3*17  arid   integrate 
to  determine  the  expression  for  flow  rate. 


Q'  -  -  ih[  £* x3  - cX  3.21 


Rearrange  and  solve  for  the  pressure  gradient. 


^T  -  ~  -L2#(Q'  +  cX)  3.22 

da       x3 


Now  solve  for  the  time  average  flow  rate  in  the  fixed  frame 
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with  the  ultimate  object  of  obtaining  an  expression  for  the  flow 
rate  in  the  moving  frame  for  substitution  in  equation  3.22. 


<Q>  =  <  I  v^  dx  > 


=  <  l   (v'  +  c)  dx  y 


3.23 


3.24 


The  first  part  of  the  integral  is  Q'.  The  second  is  the 
phase  velocity  times  R,  the  time  average  value  of  X. 


<Q>    =  Q'   +  cR 


3.25 


Rearrange. 


Q*    =   <Q>    -  cR 


3.26 


which  is  independent  of  longitudinal  position.   Substitute  Q* 
in  equation  3.22  and  note  that  c  (X  -  R)  =  cr. 


dP 
dz' 


12  ^C 
v3 


(<Q>  +  cr) 


3.27 


The  expression  for  pressure  from  equation  3.10  can  be  dif- 
ferentiated with  respect  to  z,    since  -r-  =  j~r  »  "to  yield 


12M  (<Q> 
X^ 


cr) 


d_ 
=  dz 
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f 
e 

wl 


K  (Xj  -  X) 

wl 


3.28 


For  a  given  periodic  driving  current,  equation  3.28  represents 
a  non-linear  equation  which  can  be  solved  for  X.   Only  one  boundary 
condition  is  required.   The  response  desired  is  periodic  with  the  ' 
same  wavelength  as  the  driving  current.  The  required  solution 
for  X  is  found  by  adjusting  the  boundary  condition  at  one  position 
such  that  the  deflection  returns  to  the  same  value  after  one 
wavelength,  that  is  ,  X(z')  =  X(z'  +  A). 

3.3  Determination  Of  Pressure-Flow  Characteristics  For  A  Known 
Transducer  Actuator  Deflection 

3.3 »1  Description  Of  Peristaltic  Traveling  Wave  Deflections 

The  peristaltic  pump  devel.oped  for  this  thesis,  as  stated 
in  section  3 • 1 «3 »  operates  so  that  electric  forces  dominate.   The 
traveling  wave  frequency  is  the  same  as  the  applied  electric  exciter 
frequency.   The  transducers  come  very  close  to  operating  as  open- 
close  relays.  Thus,  an  approximate  description  of  the  tube  deflection 
is  known.   The  non-linear  equation  3«28  need  not  be  solved  to 
determine  X,   Rather,  the  pressure-flow  characteristic  for  the 
pump  can  be  determined  by  integrating  equation3.27  under  the  assump- 
tion that  all  quantities  inside  the  integral  are  known  functions 
of  z. 

Figure  J>.h   is  a  Sanborn  Recording  of  the  tube  deflections 
during  pumping.   This  figure  can  be  compared  to  Figure  2.15c,  a 
similar  trace  except  that  no  fluid  filled  tube  is  in  the  transducers. 
The  erratic  traces  on  Figure  3A   are  due  to  sensitivity  adjustment 
difficulties  on  the  recorder,  not  fluid  action.  The  presence  of 
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SANBORN  RECORDER  PLOTS  OF  TRANSDUCER  ACTUATOR  MOTION  VERSUS  TIME 

DURING  PERISTALTIC  PUMPING 
Idc  B  10.0  @  26  turns    Frequency  =0.33   Ho  =  10  cm. 
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fluid  does  not  significantly  change  the  traveling  wave. 


Ap 


A  AV 


r^dp 

Jo  dz 


dz 


-  12 -a 


X 


*    <:q> 


Q>+  cr 


dz 


3.29 


3.30 


The  peristaltic  traveling  wave  can  be  approximated  as  a 
square  wave  or  as  a  sine  wave.  Both  approximations  are  made 
in  succeeding  sections, 

3.3.2  Square  Traveling  VJave  Approximation  To 
Peristaltic  Pumping 

The  peristaltic  traveling  wave  can  be  approximated  as  a 
square  wave.  This  is  illustrated  in  the  laboratory  and  wave 
frames  in  Figure  3 .5. 

The  pressure  rise  over  a  wavelength  for  a  square  wave  is 
determined  geometrically  using  equation  3 .30. 


Ap^  =  -  124 


<Q>  +  cr    <Q>  -  cr 

+ 


(R  +  r)3     (R  -  r)3 


2 


3.31 


Normalize  this  equation  by  letting 


h  = 


R 


3.32 


<Q>    - 


<Q> 
Re 


Ap 


R     APx 

GAJLkc 
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3.33 
3.3^ 


SQUARE  WAVE  PERISTALTIC  PUMPING 


x'  i 


a.  Square  wave  deflection  in  the  moving  frame.  The  wall  velocity 
is  -c  and  the  volume  flow  rate  is  independent  of  z', 


x  A 


b.  Square  wave  deflection  in  the  fixed  frame.  The  increasing 
volume  in  the  half  wavelength  between  ife   and  \     necessitates 
a  net  flow  into  that  volume.  The  velocity  profiles  would  yield 
a  net  flow  out  of  the  decreasing  volume  between  0  and  \/2. 


Figure  3*5 
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Substitute  these  values  in  equation  3.31. 


-Ap  - 


<Q>  4-  h 
(1  +  h)3- 


<Q>  .-  h 
(1  -  h)3 


3.35 


The  normalized  pressure  rise  as  a  function  of  normalized 
time  average  flow  rate  is  plotted  in  Figures  3.6,  3.7,  3.8,  and 
3.9 -for  various  values  of  h.   As  h  approaches  one,  the  channel 
is  cutoff  and  the  pressure  rise  becomes  infinite.   Such  a  positive 
displacement  pump,  operating  so  the  channel  is  cut  off  while 
pumping,  can  theoretically  pump  against  an  infinite  head. 

When  h  is  zero,  there  is  no  channel  wall  deflection.   If 
a  floitf  rate  existed  in  the  channel,  this  case  would  give  the 
pressure  rise  up  the  channel  for  a  given  flow  rate  with  no  wall 
deflection.  However,  in  the  specific  case  of  the  peristaltic 
pump  where  all  flow  is  due  to  wall  motion,  the  flow  rate  and 
pressure  rise  would  also  be  zero. 

The  intercept  points  for  the  normalized  pressure-flow  curves 
are  determined  by  evaluating  the  normalized  pressure  rise  for 
zero  flow  rate  and  the  normalized  time  average  flow  rate  for 
zero  pressure  rise.   Intercept  values  for  operating  suction  heads 
and  normalized  channel  x^idth  variations  are  shown  in  Table  3»2. 


AP, 


<Q> 


1 


1 


(1  -  h)3     (1  +  h)3 

(1  -f  h)3  -  (1  -  h)3 
(1  +  h)3  +  (1  -  h)3 
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3.36 


3.37 


Table  3.2 
Square  Wave  Pressure-?3_ow  Intercepts 


Figure 

H 
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h 
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<Q>0 

3.7 
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.049 

3.8 

20 

.172 

.196 

.082 

3.9 

30 

.214 

.321 

.122 

3.10 

4o 

.228 

.374 

.137 

3.3.3  Sinusoidal  Traveling  Wave  Approximation  To 

(17) 
Peristaltic  Pumping 

The  peristaltic  traveling  wave  can  also  be  approximated  as 

sinusoidal.  This  is  illustrated  in  Figure  3.6.  The  pressure 

rise  over  a  wavelength  for  a  sinusoidal  traveling  wave  is 

determined  by  integration  of  equation  3.30. 


Ap 


12  JUL     (^   <Q>  +  cr  cos(wt-kz)  ,     .     .   *   0 
-y—  /  ~ ~      -* zi     (  -k  dz)  3. 

K  Jo    [R  + 


r  cos  (wt  -  kz)J 


38 


,(18) 


When  solved    ,  this  integral  becomes 


A V  tf^ylZ     [3cRr2-<Q>(2a2+r2) 


] 


3.39 


Normalize  equation  3«39  using  the  same  normalization  as  in 
the  previous  section.   Then 


r 


Ap 


(1  -  h*) 


275/2 


3  - 


2_j_hJ 
,  2 


<Q> 


3.40 
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SINUSOIDAL  WAVE  PERISTALTIC  PULPING 


a.  Sinusoidal  wave  deflection  in  the  moving  frame.  The  wall 
velocity  is  -c  and  the  volume  flow  rate  is  independent  of  z'< 


/a 

b.  Sinusoidal  wave  deflection  in  uhe  fixed  frame.  The  increasing 
volume  in  the  half  wavelength  between  \/2  and  X   necessitates 
a  net  flow  into  that  volume.  The  velocity  profiles  would  yield 
a  net  flow  out  of  the  decreasing  volume  between  0  and  \/2. 


Figure  3.6 
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The  discussion  in  the  previous  section  is  applicable  to 
the  sinusoidal  approximation  of  the  peristaltic  traveling  wave. 
As  from  the  previous  section,  the  intercept  points  for  Figures 
3.7»  3.8,  3.9.  and  3.10  are  below.  Specific  values  of  the  points 
are  in  Table  3.3. 


APr 


<Q>0  = 


(1   -hV72 
2 


2  +  ti 


Table  3.3 
Sinusoidal  Wave  Pressure-Flow  Intercepts 


3. 41 
3.42 


Figure 
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3.3.4  Experimental  Pump  Operation 
The  pump  was  operated  peristaltically  with  the  same  suction 
head  and  normalized  channel  width  combinations  for  which  pressure- 
flow  characteristics  were  determined  in  sections  3»3«2  and  3.3*3 • 
Values  of  pump  flow  rates  and  pressure  heads  were  measured  and 
are  recorded  in  Appendix  I.  The  operating  points  are  pl.otted 
on  the  appropriate  pressure-flow  curve,  Figure  3.7»  3.8,  3.9>  or  3«10. 

74. 


NORMALIZED  PRESSURE-FLOW  CHARACTERISTIC  FOR  PERISTALTIC  PUMPING 

h  -  0.131 

H  =  0.1  m. 
s 


0.10 


0.05  -, 


Square  Wave  Characteristic 


Sinusoidal  Wave  Characteristic 


Figure  3,7 
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NGj'jlALIZED  PRESSURF-FLCW  'CHARACTERISTIC.  FOR  PERISTALTIC  PUMPING 
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Figure  3.8 
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NORMALIZED  PRESSURE -FLOW  CHARACTERISTIC  FOR  PERISTALTIC  PUMPING 
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Figure  3,9 
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NORMALIZED  PRESSURE-FLOtf  CHARACTERISTIC  FOR  PERISTALTIC  PUMPING 
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Figure  3t'10 
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In  Appendix  I,  the  measured  flow  rpte  was  divided  by  the 
channel  width  prior  to  normalization  to  obtain  flow  per  unit 
channel  width.  The  pump  pressure  differential  was  determined 
by  measuring  a  static  height  difference  between  the  pump  suction 
and  discharge  headers.  The  normalized  flow  and  pressure  differ- 
ential are  determined  using  the  same  definitions  as  in  sections 
3*3*2  and  3*3*3  and  applying  the  measured   parameters  of  frequency, 
wavelength,  channel  width,  and  phase  velocity. 
3.3.5  Cut  Off  Or  Occluded  Pumping 

As  vias   pointed  out  in  section  3*3*2,  when  h  =  1,  the  channel 
is  cut  off  during  pumping  and  theoretically  the  pump  can  pump 
against  an  infinite  head.  In  this  case,  a  pressure-flow  char- 
acteristic is  meaningless. 

Theoretically,  the. entire  contents  of  the  pump  are  pushed 
forward  at  the  phase  velocity  of  the  applied  traveling  wave  during 
occluded  pumping.  For  a  unit  tube  width, 

<Q>  =  Re  3*44 

Using  the  same  definition  for  normalized  time  average  flow 
rate  as  in  section  3*3*2, 


<Q>   -  4g>  3.45 


<Q> 

measured  3.^4-6 


" '  theoretical 
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OO  can  range  from  zero,  no  flow,  to  one,  perfect  occlusion. 
The  experimental  normalized  time  average  flow  rate  in  Appendix  II 
ranged  from  0.26  to  0,^5.   Considering  the  sources  of  error, 
which  will  be  mentioned  in  section  3.^»  the  most  important  to 
inefficient  occluded  pumping  are  improper  time  phasing  between 
actuators  and  the  fact  that  the  pump  consists  of  only  one  wavelength. 
Both  of  these  items  bring  about  physical  conditions  under  which 
the  entire  contents  of  the  tube  under  each  actuator  are  not  pushed 
f orward  as  the  actuator  squeezes  downward. 

If  the  pump  i^ere  operated  so  that  each  transducer  was  on  or 
off,  almost  ideal  occlusion  could  be  obtained  if  each  actuator 
was  held  down  until  the  succeeding  one  was  closed.  Under  these 
circumstances,  the  fluid  has  no  choice  except  to  move  onward. 

For  the  specific  pump  in  this  thesis  and  under  the  type 
of  operation  finally  used  where  maximum  current  was  applied  for 
all  pumping, peristaltic  and  occluded,  occlusion  is  the  most 
efficient  form  of  pumping.  For  a  given  power  input,  occlusion 
gives  the  greatest  flow  rate.  Unless  less  than  fully  occluded 
pumping  is  desired,  the  most  effective  pump  would  be  built  to 
operate  fully  occluded. 

3«3«A  Reynold's  Number 

The  appropriate  Reynold's  number  describing  flow  in  the 

(19) 

peristaltic  pump  has  been  shown  to  be 


\r  A  ■ 
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The  Reynold's  number  varies  between  1.31  and  8.65  for  the 
operating  points  of  the  pump.  Exact  valuos  are  shown  in  Appendices 
I  and  II.   In  each  case,  inertial  forces  are  greater  than  viscous 
forces.  In  spite  of  this,  the  viscous  model  developed  yields 
reasonably  accurate  results, 
3.^-  Sources  of  Error 

It  is  difficult  to  attribute  the  error  between  theoretical 
pressure-flow  characteristics  and  experimental  operating  points 
to  any  single  cause.  The  inertial  nature  of  the  pump,  whi ch  is 
not  accounted  for  in  the  theoretical  model ,  could  cause  some 
error.  The  mechanical  aspects  of  pumping  could  have  produced 
errors.  The  rubber  pump  tube  was  somewhat  elastic  which  could 
.  have  altered  flow.  The  tube  was  eliptical  in  shape  rather  than 
a  square  channel  as  the  theory  assumes.  The  model  assumes  the 
vertical  \-jalls  of  the  tube  are  always  vertical  and  perfectly 
elastic  to  up  and  down  motions,  while  the  channel  is  of  constant 
width.  The  real  pump  has  elastic  vertical  walls  which  bow  out, 
Channel  width  varies  as  the  pump  actuator  moves  up  and  down.  As 
was  pointed  out  in  section  2.5»  perfect  ninety  degree  phase  shifting 
between  transducers  was  not  obtained.   The  uneven  actuator  motion 
could  have  affected  flow.  The  pump  pressure  head,  from  which 
A?  results,  is  on  the  order  of  one  millimeter.  Accurate  mea- 
surement of  these  differentials  was  difficult.  The  pump  has 
rather  long  sections  of  flexible  tube  on  the  inlet  and  outlet. 

The  flow  resistance  and  surface  tension  effects  at  the  ends  of 
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these  tubes  could  alter  pressure-flow  relations. 

Considering  the  sources  of  error,  it  is  surprising  that  the 
operating  points  are  as  close  to  the  theoretical  press\ire-flow 
curves  as  they  are.   But,  the  real  value  of  the  project  is  not 
the  relative  agreement  between  theory  and  experiment.  Its  value 
is  in  the  physical  realization  that  longitudinal  peristaltic 
type  flow  can  be  produced  by  purely  vertical  transducer  motion. 
No  longitudinal  mechanical  motion  is  involved. 
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CHAPTER  IV 
CONCLUSIONS 

The  experimental  work  done  in  this  thesis  has  shoi^n  that 
the  observed  peristaltic  pump  operating  points  lie  reasonably 
close  to  the  theoretical  pressure-flow  characteristics  of  the 
pump.   It  has  not  been  shown  that  the  pressure-flow  characteristics 
accurately  predict  pump  operation.  First,  very  few  operating 
points  were  measured.  No  real  trend  or  correlation  between 
measured  and  theoretical  operation  could  be  observed.  Second, 
the  large  number  of  error  sources  make  it  difficult  to  differentiate 
flow  variations  due  to  theoretical  prediction  from  flow  variations 
due  to  experimental  error.  The  exception  to  this  is  occluded 
pumping.   Complete  occlusion  theory  predicts  a  normalized  flow 
rate  of  one.  Flow  rates  less  than  one  are  due  to  inefficient 
or  improper  experimental  operation  of  the  pump  and  are  readily 
recognizable. 

The  thesis  and  experimental  work  have  been  of  value.  The 
major  reason  for  the  effort  was  to  show  that  peristaltic  pumping 
could  be  produced  by  a  principle  which  had  never  been  tried. 
Longitudinal  flow  was  produced  where  no  longitudinal  mechanical 
motion  was  involved.  The  importance  of  this  fact  cannot  be  over 
emphasized.   The  longitudinally  traveling  peristaltic  wave  was 
produced  by  properly  phasing  the  transverse  motion  of  discrete 
electromechanical  transducers.  If  a  peristaltic  pump  is  ever 
to  be  miniaturized  and  used  as  a  replacement  for  an  organ  in  the 
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body,  the  use  of  some  sort  of  mechanical  roller  is  not  practical. 
But,  such  a  pump  could  be  designed,  using  the  principles  of  this 
thesis,  which  would  be  easy  to  miniaturize,  mechanically  simple, 
flexible,  light  weight,  compact,  and  reliable.  For  instance, 
if  very  small  displacements  are  desired,  the  transverse  transducer 
motion  could  be  produced  by  piezoelectric  expansion  of  ceramic 
crystals.   The  excitation  for  a  piezoelectric  crystal  is  a  voltage 
applied  across  it. 

Another  new  principle  developed  in  this  work  is  the  direct 
relationship  between  the  applied  electric  force  and  the  transverse 
motion  of  the  transducer  actuator.   In  the  case  of  the  peristaltic 
pump  developed  in  chapters  two  and  three,  the  magnitude  of  the 
electric  force  is  inversely  proportional  to  the  gap  width  squared. 
In  future  peristaltic  devices,  it  might  be  desirable  to  utilize 
the  electric  force  or  the  magnetic  field  directly  to  produce 
fluid  motion.   This  has  possibilities  in  pumping  magnetic  liquids 
directly  or  in  pumping  by  means  of  magnetic  liquids.  A  traveling 
magnetic  wave  might  be  made  to  pull  magnetic  liquid  along  just 
as  the  traveling  field  of  a  motor  pulls  the  armature  along.   Both 
of  these  principles  are  of  value  for  future  experimental  work. 

The  specific  pump  developed  could  be  improved  in  the  following 
ways.  For  most  efficient  operation,  it  should  pump  by  totally 
cutting  off  the  flow  channel  in  its  transverse  motions.   This 
causes  the  greatest  volume  flow  rate  for  power  input.  However, 
total  occlusion  could  be  disadvantageous.   Blood  is  very  sensitive 


to  pumping.  The  effect  of  being  pinched  in  the  flow  channel  could 
damage  blood  so  that  it  might  not  be  usable.   The  pump  should 
be  improved  for  pumping  both  fully  occluded  and  peristaltically. 
Phasing  and  phase  sequencing  need  improvement.  Proper  application 
of  excitation  could  produce  rhythmic  heart-beat-like  pumping 
rather  than  continuous  smooth  flow,  if  that  were  desired.  Since 
the  transducers  act  almost  like  on-off  relays,  a  complex  sinusoidal 
excitation  is  not  needed.  More  efficient,  simpler  excitation, 
which  would  be  just  as  effective,  could  be  provided  by  a  DC  on-off 
power  source. 

One  of  the  most  serious  defects  in  the  pump  is  the  low  volume 
flow  rate.  The  pump  tube  used  was  a  long,  thin  balloon.  It  vras 
difficult  to  find  an  elastic  tube  which  would  fit  in  the  gap 
between  the  actuator  and  the  exciter  and  which  did  not  exhibit 
elastic  forces  opposing  gap  close  off  greater  than  the  electric 
forces  applied.   Ideally,  a  long,  very  flexible,  broad  tube  is  needed 
which  would  fill  more  of  the  channel  width  and  have  negligible 
vertical  wall  rigidity.  Total  channel  "elasticity"  should  be  due 
to  the  suction  height  of  the  fluid  in  the  channel. 

The  type  of  electromechanical  peristaltic  pump  investigated 
in  this  thesis  is  rugged,  inexpensive,  and  simple.  It  was  built 
using  common  transformers,  springs,  plastic,  and  balloons.  The 
pump  operated  consistently  and  dependably.  These  features  ,  and 
those  already  mentioned,  make  it  worth  additional  development. 
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